New class of level statistics in correlated disordered chains 
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We study the properties of the level statistics of ID disordered systems with long-range spatial 
correlations. We find a threshold value in the degree of correlations below which in the limit of 
large system size the level statistics follows a Poisson distribution (as expected for ID uncorrelated 
disordered systems), and above which the level statistics is described by a new class of distribution 
functions. At the threshold, we find that with increasing system size the standard deviation of 
the function describing the level statistics converges to the standard deviation of the Poissonian 
distribution as a power law. Above the threshold we find that the level statistics is characterized 
by different functional forms for different degrees of correlations. 
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Investigating statistical properties of the energy spec- 
trum and the behavior of the level statistics proved to 
be a useful approach to study electronic properties of 
disordered systems 0, 0, In 3D disordered systems 
the Anderson transition between metallic and insulating 
phase is associated with a transition in the level statis- 
tics distribution from a Wigner-Dyson form to Poisson 
distribution 0, [H El j an d the level statistics at the criti- 
cal point has been studiedQ. In ID disordered systems, 
for any degree of disorder, the level statistics is described 
by the Poisson distribution in the limit of large system 
size corresponding to electronic localization and insulat- 
ing behavior 0- 

Recently, it has been demonstrated numerically that 
introducing long-range correlations in the spatial order of 
atoms with different energies in a chain can lead to elec- 
tronic derealization creating an interesting debate^]- 
Further, it has been shown that there is a localization- 
delocalization transition at a critical value of the degree 
of correlations imposed on the disorder in the system, 
and later works extended these results to other mod- 
els 0, 0, 0, 0| . Also, this type of transition can 
be found in quasiperiodic systems, as the Aubrey- Andre 
model and other models 14'], indicating the importance of 
some type of ordering. These theoretical findings are sup- 
ported also by experimental results showing derealiza- 
tion and electronic transport driven by extended states in 
correlated-disordered GaAs/Gao.7Al .3 superlattices (15) . 
Recently, a metal-insulator transition has been reported 
in 2D correlated disordered systems [lqj. 

Here we hypothesize that, as the localization proper- 
ties of the electronic states in a disordered system are 
affected by the degree of spatial correlations || El, the 
properties of the level statistics of the energy spectrum 
of such correlated disordered systems could also change. 
The effect of correlations in the level statistics can be 
combined with the fact that interactions, which are ex- 
tremely important in ID, also affect the level statistics 
[13. Specifically, we investigate how the functional form 
of the distribution describing the level statics is affected 
by the degree of correlations introduced in the system. 



We demonstrate that the Poissonian form describing the 
level statistics in disordered ID systems in the thermo- 
dynamic limit is preserved even when certain degree of 
spatial correlations is introduced. Further, we find a criti- 
cal threshold for the degree of correlations above which in 
the thermodynamic limit there is a transition to a differ- 
ent class of distribution functions for the level statistics. 

We consider the standard ID tight-binding hamilto- 
nian with nearest-neighbor interaction 
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where V is the coupling energy and i ranges from 1 to 
N, where N is the system size. To fix the energy scale 
we choose V = 1. In the case of uncorrelated disorder, 
the site energies are randomly drawn from a certain 
probability distribution, commonly a box (uniform) dis- 
tribution or a Gaussian. This is equivalent to consider 
the series of site energies as white noise. In contrast, 
for systems with correlated disorder we introduce spa- 
tial long-range correlations in the series of site energies 
so that their sequence describes the trace of a frac- 
tional Brownian motion. To this end, we obtain the site 
energies using the inverse Fourier transform 



N/2 

E 

k=l 



.-0 



1-/3 



1/2 



N 



(2) 



where fa are N/2 randomphases uniformly distributed 
in the interval [0, 2tt] Thus by construction the 

power spectrum of the series {&} is of the type l/fc' 5 . By 
choosing different values for the exponent we generate 
series of site energies with different degree of spatial cor- 
relations: for (3 = we have pure disorder (white noise), 
[3 < corresponds to anticorrelations, and [3 > rep- 
resents positive correlations in the series of site energies 
In our study we consider only systems with positive 
correlations ((3 > 0). 

Once the series {^} is obtained, we normalize it to zero 
mean and unit standard deviation, thus fixing the width 
of the site energy distribution to unity. This is equivalent 
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to keeping the 'traditional disorder' of the system fixed, 
since the standard deviation of the distribution of 
quantifies the variety in the site energies of the atoms 
forming the chain, while their spatial order is quantified 
by the exponent j3. Note that for uncorrelated-disordered 
systems (/3 = 0) the disorder is quantified by the standard 
deviation in the case of site energies randomly drawn 
from a Gaussian distribution, or by the width of the box 
in the case of a box distribution. 

After the normalization of the site energies, we diago- 
nalizc the hamiltonian flf to obtain the energy spectrum 
{Ei}, where E\ < E% < . . . < E^. For any system size 
N and any value of the correlation exponent (3 in our nu- 
merical calculations, we diagonalize 2 2i /N realizations 
of the hamiltonian (|TJ . Thus we have a sufficiently large 
ensemble of realizations to avoid statistical fluctuations 
in our results, while we consider the same number of 2 24 
energy levels for any N. 

Once the energy spectrum is obtained, we study the 
distribution of the spacings between consecutive energy 
levels. Since the density of energy levels is not constant 
throughout the energy band, and thus the local average 
energy spacing is not constant either, one cannot com- 
pare fluctuations in the spacings obtained from different 
regions of the band. To avoid this problem, we normalize 
to unity the local average energy spacing from different 
regions of the energy band, thus effectively normalizing 
all energy spacings to the same scale. This 'unfolding' 
of the energy spectrum is a procedure commonly used 
in the study of level statistics of disordered systems [2(| . 
In brief the unfolding procedure consists of the following 
steps: we first introduce the integrated density of energy 
levels g{E) defined as 



l(Ei) = i. 



(3) 



Thus, g{Ei) is the number of energy levels below the en- 
ergy Ei. Second, we fit g(E) using a polynomial function. 
This fit represents the averaged integrated density of en- 
ergy levels g{E). Next, the unfolded energy spectrum 
{ei} is obtained from the map 



g{Ei). 



(4) 



To avoid unfolding problems related to irregular be- 
havior of g(E) at the borders of the energy band, we 
consider only energy levels from the central region of the 
band. Specifically for a system of size N, we obtain N 
energy levels (eigenvalues of the hamiltonian (1)), and 
we consider the central part of the spectrum {Ei}, where 
i G [N/3 + 1, 2N/3]. We obtain the averaged integrated 
density of levels g(E) by fitting g(E) with a cubic poly- 
nomial in the interval [E N / 3+ll E 2 n/3\- 

Using the unfolded spectrum {ei}, we study the nor- 
malized distribution function of energy spacings P{s), 
where Si = £j+i — £j. From Eqs. 0, Q, we have that 
the average level spacing is (s) = 1. 

For the classical cases of ID and 2D uncorrelated dis- 
ordered systems in the limit of large system size P{s) 
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FIG. 1: Log-log plot of a as a function of the system size 
N for different values of the correlation exponent (3. Inset: 
same dependence in linear scales. Solid lines represent fits 
with Eqs. ©, and ©. 



follows the Poisson distribution 



Pp{s) = e 



(5) 



A Poisson distribution for P(s) indicates strong cluster- 
ing between energy levels because it reaches maximum 
when s — > (Fig. 0J. 

In our analysis we characterize P(s) using its standard 
deviation a. For the classical case of Poissonian form 
for P(s) J5}, we have a = op = 1. For convenience, we 
study the behavior of a = 1 — <r. In Fig. ^ we show a 
as a function of the system size N for different values 
of the correlation exponent (3. For uncorrelated disorder 
((3 = 0) and for large N, a — > (or equivalently a — > 1), 
indicating Poissonian behavior for P(s) as expected for 
ID disordered systems 0. 

Introducing certain degree of spatial correlations ((3 > 
0) in the system we find that in the limit of large N the 
level statistics exhibits again Poissonian behavior (p — > 0, 
Fig.^i), indicating energy level clustering and electronic 
localization. This finding indicates that Poissonian be- 
havior in the level statistics exists even in the presence of 
long-range correlations. However, we find that the con- 
vergence of P{s) to the Poisson distribution with increas- 
ing N is much slower in the case of correlated disorder 
((3 > 0) compared to the case of uncorrelated disorder 
(P = 0)(Figs. [HandEJ. 

Increasing the strength of the long-range correlations 
we find a threshold value, a critical exponent (3 C , above 
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which the functional form of P(s) does not converge to 
Poisson distribution — i.e., b does not converge to zero 
in the limit of large N (Figs. and|2J). At the critical 
value C we find that for increasing N the level statistics 
converges very slowly to Poisson distribution and b — > 
as a power law (Fig.QJ)). 

We next investigate the functional dependence of b on 
the system size N for the three regimes: (i) < C ; (ii) 
— C and (iii) > C . Based on simulations of various 
systems sizes up to N — 2 17 we model the behavior of b 
using the following expressions: 



ai7v -6 1 -c 1 logAr ^ /3< c 

a 2 N- b2 (0 = C ) 

boo + a 3 N~ b3 [ft > C ) 



(6) 
(7) 
(8) 



where all the parameters are positive and in general de- 
pend on the correlation exponent 0. To test the validity 
of these expressions, we use the Levenberg-Mardquardt 
(L-M) algorithm (|2j) to fit the data in Fig. [T]a,b and 
to estimate the optimal values of the parameters. Al- 
though we have tried more expressions to fit the data, 
those in Eqs. JO}, Q and © are the best we have found 
attending two criteria: good description of data, and less 
number of parameters. 

We note that for (3 < f3 c the behavior of b vs. N on 
log-log plot presents negative curvature (Fig. QJi), while 
for (3 — f3 c , b depends on TV as a power-law with nega- 
tive slope, so that in both cases for N — > oo, b — > 
(Fig. \I]p) indicating that in the thermodynamic limit the 
level statistics is Poissonian. In contrast, for (3 > (3 C , the 
dependence of b on N has positive curvature (FigHJ)), 
indicating a decay slower than a power-law. We find 
that data are best modeled (with best fit based on the 
L-M algorithm) as a power law with an additive positive 
constant boo such that b = boo for N — > oo (JSJ. Thus 
our results suggest that for (3 > C , in the thermody- 
namic limit b does not converge to zero, and that the 
level statistics is not of Poissonian type. 

To determine f3 c we use the following procedure: (i) 
starting from small values of (3 for which the dependence 
of b on N follows Eq. @), we increase and observe that 
the fitting parameter c\ decreases, and for a given value 
of (3 becomes zero, so that Eq. © is not valid anymore; 
(ii) starting from large values of (3 for which b follows Eq. 
© , we decrease (3 and observe that the fitting parameter 
boo decreases, and for a given value of becomes zero, 
so that Eq. (JSJ is not valid anymore. We find that for 
both (i) and (ii) the transitions occur at a critical value of 
/3 = f3 c = 1.55 ± 0.05 (Figs. [I]and|2J, where the behavior 
of b is described by Eq. J7J . 

We next obtain a phase diagram of the properties of 
the level statistics as a function of the degree of spatial 
correlations in the system. We systematically investigate 
the asymptotic behavior of b in the limit of large system 
size N as a function of the correlation exponent (3 (Fig. 
|2J) . For a fixed system size N we calculate how b depends 
on the spatial correlations choosing a dense set of (3 val- 
ues. We then repeat the calculations for increasing N. 
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FIG. 2: Dependence of a on the correlation exponent j3 for 
varied system size N. The solid line represents the behavior of 
<7 in the limit of TV — > oo (i.e., a x ) — a phase transition from 
Poissonian (a = 0) to non-Poissonian (a 7^ 0) level statistics 
at the critical value f3 c . 



We find that b is an increasing function of (3, and that 
for each value of N a relatively flat region at small (3 is 
followed by a sharp increasing in b for large (3 (Fig. 0). 
This behavior becomes more pronounced with increasing 
N. Further, we find that the flat region in b extends 
to intermediate values of (3 and rapidly approaches zero 
with increasing N. This is in agreement with our finding 
of level statistics of Poissonian type even in the presence 
of a moderate degree of spatial correlations in the system, 
and with the predictions of Eqs. ^ n contrast, for 

large values of (3, the values of b remain large and do not 
decrease substantially with increasing the system size N . 
Thus, we observe a transition in b centered at intermedi- 
ate values of 0, which becomes more abrupt with increas- 
ing N. To extrapolate the behavior of the level statistics 
in the thermodynamic limit and for large (3, we use Eq. 
© since for N — » 00, b — > boo- We estimate boo for a 
dense set of large and decreasing values of (3 — the solid 
thick line in Fig. [3 which sharply decreases to boo = 
for (3 = 0c = 1.55. This suggests a phase transition from 
a Poissonian behavior of the level statics characterized 
by CToo = for j3 < j3 c , indicating strong clustering be- 
tween energy levels, to a non-Poissonian phase defined by 
boo ^ for (3 > (3 C . As boo is a function of (3, this sug- 
gests that for any (3 > (3 C a different level statistics P(s) 
is obtained. Thus we find a new class of correlated dis- 
ordered systems characterized by energy level repulsion, 
different values of boo and different distribution functions 
for the energy spacings P(s). 

Similar conclusions can be drawn with the study of 
S 2 (L), i.e., the variance of the number of levels in boxes 
of length L in the unfolded spectrum. It is known that 
for Poissonian behavior, S 2 (L) is linear with L with slope 
1. We obtain this linear behavior for any j3 < (3 C (see the 
case = 1 in FigEJ), indicating Poissonian behavior, in 
agreement with our previous results. For > C , we 
obtain a non- linear and slow increasing of S 2 (L) as a 
function of L (slower for increasing 0), indicating level 
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FIG. 3: Behavior of E 2 as a function of L for spectra obtained 
for different values of and for a system size N = 2 17 . 
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FIG. 4: Distribution functions P(s) obtained for different val- 
ues of and for N — 2 17 . Inset: log-linear plot of P(s) for 
two values of and for Poisson and Wigner-Dyson (WD) dis- 
tributions. For = 1 < C P( s ) is Poissonian, while for 
= 2.4 > C the tail of P(s) decays even faster than for WD. 



repulsion and non-Poissonian behavior, also in agreement 
with our previous results, and with the behavior of the 
P(s) functions (see below). 



In Fig. 2] we show P(s) for several values of and for 
finite system size N = 2 16 . We note that in the thermo- 
dynamic limit of large N the form of P(s) may change. 
However, this change is expected to be not significant, 



since for N = 2 lb , a 



0.005. In general, we 



obtain that for < C , although there exists a moderate 
degree of correlations in the system, P(s) is exponen- 
tial (Poissonian) (see the case = 1 in FigQJ. When 
> C , the Poissonian behavior is lost: as departs 
from C , the functions P(s) for low s decrease gradually, 
and simultaneuously, an increasing peak for increasing 
appears at s = (s) =1, indicating strong level re- 
pulsion. For extreme values of /3, P(s) for low s is very 
small, while the peak at s = 1 becomes huge (see the 
case of = 3 in FigQJ. This finding is consistent with 
the expectation that extreme values of correspond to 
an ordered system, for which the level statistics is of the 
type P(s) = 8(s - 1). 

In summary we find that introducing spatial long- 
range correlations in ID disordered systems leads to a 
transition from a Poissonian to a new class of functional 
forms describing the level statistics in the thermody- 
namic limit. Further, we find a critical value for the 
correlations below which the level statistics exhibit Pois- 
sonian behavior associated with energy level clustering, 
similar to the one observed in uncorrelated disordered 
systems. Above this critical value the system is charac- 
terized by level repulsion. These findings may relate to 
previous reports on localization-delocalization transition 
in the electronic properties of ID systems driven by spa- 
tial correlations in the disorder ||. In that work, the 
transition is detected at = 2 > C , where we already 
observe non-Poissonian behavior, as expected in the ex- 
tended regime. 
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